Journal of Research of the National Bureau of Standards Vol. 60, No. 2, February 1958 Research Paper 2824 

Energy Spectra of Cascade Electrons and Photons 1 

Charles A. Olson 2 and Lewis V. Spencer 

The equations for the energy spectra of electrons and photons in a cascade shower are 
written in a form suitable for numerical applications with accurate cross sections. Trial 
calculations were carried out to check the feasibility of a step-by-step numerical integration 
procedure similar to that used successfully at low energies in noncascade problems. Numerical 
results were obtained for lead by using a monoenergetic source of electrons at 360 mc 2 and 
extreme relativistic cross-sectional forms, which assume complete screening. The method 
was designed to permit the determination of cascade-shower spectra, spatial distributions, 
and directional distributions that are exact in the sense that the main limiting factor is the 
accuracy with which the cross sections are known. The numerical results of the trial calcula- 
tions were found to agree with previously published results based on more restrictive methods 
of analysis. 

1. Introduction 

High-energy electrons and photons dissipate their energy by generating secondary electrons 
and photons. The complex radiation arising in this manner from a high-energy electron or 
photon source is called a cascade shower. Most attempts to describe this phenomenon have 
used approximations thai are valid only at very high energies, i. e., in the extreme relativistic 
range; and the results have only limited applicability in a broad range of lower energies where 
the cascade phenomenon is still of importance. At these lower energies realistic calculations 
must be largely numerical because the accurate cross-sectional forms are too complicated to 
permit analytic transformations. 

In this paper a calculations! procedure which applies without restrictions on cross-sectional 
forms, and which has been used in noncascade problems, is adapted to the cascade situation. 
It should provide an adequate way of describing cascade-shower spectra in any energy range. 
Because the equations are part of a linked set of moment equations, which describe also spatial 
and directional distributions, and possibly fluctuations, the procedure in this paper may 
eventually form the basis for a more complete description of many aspects of the cascade 
problem. 

The equations are applied, in a trial calculation, to the special case of a monoenergetic 
electron source. Calculations for an electron source are more important than for a photon 
source for the following reason: If a photon source is assumed, it is possible to calculate the 
resulting photon spectrum at all energies. This type of calculation has been done many 
times and offers no difficulty (for an example, see reference [l]). 3 One can then determine by 
simple quadrature the electron spectrum generated by this photon spectrum. The whole 
cascade spectrum of electrons and photons, except for the initial photon component already 
determined, can then be determined by appropriate superposition of solutions for mono- 
energetic electron sources. Thus, calculations for a monoenergetic electron source demonstrate 
that the main difficulty blocking a complete numerical analysis of the cascade problems is 
removed. 

For the initial calculation, extreme relativistic cross-sectional forms that assume complete 
screening are used. They include the same types of singularity as the exact cross sections 
and therefore provide a meaningful trial run. Also they are quite similar to the cross-sectional 
forms used in earlier applications of conventional cascade theory; so they make possible a direct 
check with previously published results. 

i This work was supported by the Office of Naval Research and the Atomic Energy Commission. 

2 Submitted to the University of Maryland in partial fulfillment of the requirements for the M. S. degree. 

3 Figures in brackets indicate theli'" ture references at the end of this paper. 
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2. Interactions Involved in Cascade Showers 

The two types of interactions that predominate in the development of cascade showers 
are bremsstrahlung and pair production. In a bremsstrahlung interaction a high-energy 
electron or positron, under the influence of a nuclear or electronic field, generates a photon. 
In pair production a high-energy photon in the electromagnetic field of a nucleus or electron 
gives up all its energy to generate an electron and a positron. The combination of these two 
processes results in the creation of photons, electrons, and positrons in increasing numbers at 
lower energies. 

At all energies below the extreme relativistic range, there are additional effects that influ- 
ence the development of the showers. There are energy losses due to atomic interactions that 
slow down the electrons more or less continuously along their path. Some of these interactions 
result in secondary electrons with high energies, i. e., knock-on electrons. Compton scattering 
also contributes some high-energy electrons to the shower, but at low energies it is of primary 
importance through its influence on the photon spectrum. The photoelectric effect and 
photonuclear interactions have cross sections that are much smaller than the other interactions 
at all cascade-shower energies, though the former is of importance below shower energies. 

3. Discussion of Earlier Work 

An excellent discussion of the lines along which the mathematical description of the cascade 
phenomenon has been developed can be found in a book by B. Rossi [2]. Because the proce- 
dural methods of this paper are directed mainly toward a solution of the cascade problem at low 
energies where accurate cross-sectional forms are needed, comments here on earlier work will 
be restricted to developments pertinent to these lower energies. 

The first attempt to describe a low-energy cascade spectrum was made by Arley [3]. He 
neglected interactions that slow down the electrons, i. e., ionization losses, above the critical 
energy; 4 and he further neglected bremsstrahlung and the generation of secondary electrons 
by all processes at energies below the critical energy. His results as compared with more 
recent calculations are low by a factor of 3 or 4. 

Tamm and Belenky [4] have obtained a much more accurate solution by an analytic 
method, using asymptotic cross-sectional forms for bremsstrahlung and pair production, assum- 
ing continuous ionization loss and neglecting the Compton effect and knock-on electrons. 
Their solution has proved to be more accurate than was at first expected, considering the 
approximations. Friedman [5] has discussed their solution at length and computed a first 
correction term. This correction term never exceeded 4.5 percent of the original solution in 
the case that Friedman considered. 

Perhaps the most satisfactory treatment of the energy-spectra problem for low energies 
has been given by Richards and Nordheim [6]. They have used reasonably accurate repre- 
sentations for all important interactions, but they treat small radiation and ionization losses 
as continuous, thus incurring an error of up to 8 percent, according to their own estimates. 
They obtain results by numerical integration for energies from % to 10 times the critical energy. 
Their results are not absolute and do not describe the spectra near the initial energy. To 
obtain more complete spectra, it is necessary to normalize their equations to a solution valid 
at high energies, and they have used the results of Rossi and Greisen [7] for this purpose. 

All these calculations are difficult to generalize to a more accurate treatment of all the 
different interactions, and hence there is a need for a new approach that will permit use of 
arbitrarily accurate cross sections. 

4. Analytic Treatment of the Problem 

In this section equations will be derived for the energy spectra of electrons and photons. 
No distinction is made between positrons and electrons, the two being considered as having 

4 The critical energy in a given material is the energy at which the radiation losses of an electron equal the ionization losses. 
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identical cross sections and combining to make up an "electron" spectrum. The differentia] 
energy spectrum for electrons will be represented by 1(E) dE, which is the track length traveled 
by all electrons in the energy range E to E+dE. H(a)da represents the photon spectrum simi- 
larly. 5 In all cases energies are in units of mc 2 , and the electron energies include the rest energy. 

To set up equations in terms of 1(E) and H(a), a time-independent source is assumed. 
Physically, this means a continuous cascade in which the average number of electrons, positrons, 
and photons at any given energy remains constant with time. The equations will be derived in 
a general form, and considerations of the interactions and cross-sectional forms will be reserved 
for a separate section. 

Continuity equations for 1(E) and H(a) can be written as follows : 

r dE'k es (E' ,E)I(E')=s eg {E)+ rda'k es {a r ,E)H(a') (1) 

J E J E 

rda'k ps (a',a)H(a')=s p8 (a)+ f" dE'k pg [E' ,a)I{E') , (2) 

where 

k es (E',E) is the total probability per unit path length traveled for an electron of energy E' 
to undergo an interaction that removes it from the energy range above E, 

k eg (a,E) is the total probability per unit path length traveled for a photon of energy a 
to generate an electron of energy greater than E, 

s eg (E) is the total source strength for generating electrons with energies greater than E, and 

k ps (a',a), k pg (E',a), and s pg (a) are similar quantities with the photon and electron roles 
reversed. 

Note that the superscripts refer to "electron scattering," "photon generating," etc., al- 
though strictly speaking the interactions include types other than scattering. Physically, eq 
(1) states that for a time-dependent source the rate at which electrons are removed from the 
energy range above E is equal to the rate at which electrons come into the energy range above 
E. Equation (2) has the same meaning for photons. 

Equations (1) and (2) completely define 1(E) and H(a) for specified sources; however, 
these equations as they stand are not very suitable for numerical work. It is advantageous to 
make further transformations suggested by the nature of the variations of spectra and cross 
sections. The particular choice given here is based largely on numerical and analytic ex- 
ploratory calculations that are not sufficiently pertinent to the main discussion to warrant a 
detailed description, but some general arguments to justify this choice are given below. 

1(E) may vary rapidly if s eg (E) is discontinuous (i. e., if the source is concentrated at some 
energy), as for example, in the case of figure 1 in reference [8]. A more slowly varying function 
on which to base a step-by-step numerical integration is 



F(E)=rdE'I(E') y 



which gives the track length traveled by all electrons while their energy is above E. The func- 
tion H(a) may also be singular if the source function s vg (a) is discontinuous, but photon cross 
sections permit the removal of such singularities by iteration (see, for example, reference. [9]). 
It is therefore not necessary to utilize an integral form of the photon spectrum. 

Equation (1) could be rewritten in terms of F(E) by partial integration, but this results in 
the differentiation of k es (E / ' ,E), which is already singular as E f approaches E. A further trans- 
formation is needed to make the equation more tractable, namely integration of all terms of eq 
(1). With regard to eq (2), k ps (a,a) and k pg (E',a) are integrals of nonsingular cross sections. 
It is therefore possible, and in fact advantageous, to accompany the partial integration of the 
last term in eq (2) with a differentiation of all terms in the equation. 



6 (E) and //(a) are in the units g cm~ 2 (Mc*)-i sec -1 . These are track-length units referring to a source of one electron. The dimensions of 
1(E) and 11(a) change to those of flux if the source is 1 electron per gram per second. 
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Applying these transformations we obtain for the left side of eq (1) 

CdE' rdE'W(E",E')I(E")= r<lE'I(E') {' ' dE"V°{E',E") 

J E J E' JE JE 

= CdETiE')^ C'dE"k es {E',E"). 

The other term in the integration by parts, —F (£") dE"k e *(E'E") | , is zero at both 

LJ E A\e 



limits. 

Equation (1) thus becomes 



where 



C dE'K es {E f ,E)F{E')=S e8 (E) + f° ' da'K«(a' \E)H(a') y (3) 

J E J E 

K es (E',E)=~r dE"k es (E',E"), 
K<*{a',E)= f" dE'k es (a',E'), 

J E 

S eg (E)= fdE'^E'). 
Equation (2) becomes 

»(a)H(a)- C da'K v *{cL' ,a)H(a')=S™(a)+ C dE'F(E')K™(E' ,a), (4) 

J a J a 

Note that the other term on the right side, 

~{F(E')k^(E',a)\^ ma }, 

obtained after integration by parts, is zero at both limits. 
The various quantities in (4) are defined as follows: 

n(a)=k vs (a,a) (4a) 

is the total photon interaction probability per unit path length for photons at energy a. 

K*>(a' \a) = bk "£'" ] ] =NZa{a' », (4b) 

where a {a! ,a) is the photon differential scattering cross section per electron, and N is the 
number of atoms per gram. 

where <j(E' ,a) is the cross section per atom for producing radiation. 



*•<«> —*&■ «<« 



s the photon source strength at a. 



5. Interactions and Cross Sections 

Equations (3) and (4) as derived in section 4 make no reference to specific interactions or 
cross-sectional forms, and consequently this choice can be adapted to the specific nature of each 
problem. In terms of the interactions discussed in section 2, the functions K may be written 

xi. \lL , hi ) — xibremsstiahlung losses T~ -^-atomic interactions? 

z£ (a ,i^J =/lpair production I •^Ccmpton electrons > 

K ps (a' \ol) =i?Compton photons, 
xi. \Hj ,oO—^bremsstiahlung photons > 

M(«) === At( a )pairproduction+MWComptonphotons + M(«)absorption- 



At this point accurate cross-sectional forms could be introduced in the equations, but an 
examination of their suitability can be made with simplified forms similar to whose used by 
earlier investigators. Specifically, (a) extreme relativistic cross sections are used, (b) com- 
plete screening is assumed, and (c) the rest mass of the electron has been neglected insofar as 
it influences the limits of the integrals. Detailed cross sections are given in appendix 1. The 
K functions derived from these cross-sectional forms are 



K<\E',E)=N4${[-~ In (l->?)+! {l-r,)~ ,»-| /| In 183Z~* 
where ri=E/E'; 



where rj=E/a' ; 
where rj = a/a ; 
where rj = a/E / ', and 



8 a 
K*'(E',a)=N&±[(~ r, 2 +2rj z ) In 183Z-K-J r? 2 ], 

^a)=N^ [| In mZ-yt-^+NZ | 7p T A [ln2a+i]. 

6. Numerical Treatment of the Problem 

In determining a convenient system for the numerical integration, we note that the func- 
tions K depend primarily on ratios between the two energies involved. We wish to take 
advantage of this scaling property by choosing integration steps of equal length in the logarithm 
of this energy ratio. This choice of integration mesh has the additional advantages that one 
can integrate over an enormous energy range with comparatively few steps, and that a cross- 
sectional tabulation can be made that is applicable to widely different energy ranges. 

The numerical integration is performed by Simpson's rule, using modifications of the 
familiar (|, f , f , £, J) A set of weights to make the rule apply also to an odd number of intervals. 
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The integrals were computed essentially as follows : 

fdE'K es (E',E)F{E') = frf(ln E')E'K es {E' ,E)F{E') 

=S u ' tn E i K»(E t ,E n )F(E t ) 

1 = 

=S <* in K e *(E t ,E n )F{E { ). 

Note that the integration weights w iw decrease in size as the energy decreases. More details 
are given in appendix 2. 

Despite the transformations discussed in section 4, the function K es (E',E) in eq (3) is still 
singular. However, this singularity causes no difficulty if one writes the integral in the fol- 
lowing manner: 

rdE'K"(E',E)F(E')= [" dE'K"(E' ,E)[F{E')-F(E)]+F(E) f °dE' 'K es \E' \E) . 

Je Je Je 

Note that the first term on the right is no longer singular and can be integrated numerically. 
The last term can be evaluated without difficulty by integrating the singular parts of K es (E',E) 
analytically. Notice also that as E f ->E, the first integrand is finite, approaching the limit 
A[dF(E)lbE], where A=\im(E'—E)K e8 (E',E). This termis very small andneeds to be evalu- 

E f ->E 

atecl only approximately. A three-point numerical differentiation formula, described in appen- 
dix 2, was adopted for this purpose, the corresponding weights being denoted by a t . By 
utilizing these devices, we represent the first integral of eq (3) as the sum of terms 

F(E n ) [jT» (E n ) - § <* J5? - (Ei,En) +Aa n ] +5>*«#" (E u E n ) +Aa,] F(E t ) , 

i=0 i=0 

J 00 
dE'K es (E',E). Note that a t is zero except for i=n, n—1, n—2. 
E 

If the functional form of K es (E',E) is that given in section 5, 
jf-(£)=A^4^ ^[-^-|(l-^)ln(l-^)-H|(l-^)+^ 2 +i^]lnl83Z- 1 / 3 -^ 



~(l-v) In (i-,)+±(i--,)+^J. +NZ%p T hi "^ " 

where 7]=E/E 0) and E is the highest energy of the source of electrons and photons. For a 
discussion of the choice of the cutoff value Q, see reference [8]. In this paper it has been evaluated 
as follows: 



H(S)V-i =4 - 9214X10 ~ n ' 



where 7 =13.5 ev. 

The function K pg (E f , a) in eq (4) is not singular, but the exact expression varies rapidly as 
E'-^a. This can easily be taken into account by using addition-subtraction as before: 

J» oo r* co f* oo 

dE'K™(E',a)F(E') = dE'K™{E',a){F{E')-F{E) B=a ]+F(E) B=a \ dE'K™(E',a). 
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Written as a sum, this term becomes 

where Jf« r (a)= J dE'K* g (E' ,a) . If the functional form oiK pg (E',a) is that given in section 5, 

where rj = a/E Q . 

Equations (3) and (4) can now be written in their final form for numerical evaluation: 



1 r* n — \ ra-1 ~\ 

W=t4 -T.Wr n F{E t )+JZW\iH{cc t )+S™(E H ) 

^ \&n) L i = »=0 J 



(5) 
(6) 



where 



WVn=<*inK"{E i ,E n )+Aa t , 

Wi*=u> in K e *( ai ,E n ), 



C(E n )=jf"(E n )-j: Utn K«(E t ,E n )+Aa n> 

i=0 



(7) 



C p (a n ) = m(« n ) — co„J5>'* («„,<*«) , 



These seven functions depend only on cross-sectional forms, and they are evaluated separately 
prior to the main calculation. Note that eq (5) depends on F(E t ) and H(a t ) but not on H(a n ), 
while eq (6) depends on F(E n ) and on F(E t ) and H(a t ) . Therefore, eq (6) must be determined 
after eq (5) at each step in the integration. The problem of how to start the calculation is 
taken up in the next section. 

7. Special Solution Near E 

F{E) and H(a) must vanish at the highest energy, E . It would be possible to use eqs 
(5) and (6) , with the cross-sectional forms that have been listed, to solve for the spectra near 
E . But it was found more convenient and accurate to do a separate calculation because the 
solution near E has certain special features. In particular, we note that (1) for small energy 
losses it is always possible to use extremely simple cross-sectional forms, (2) the linkage be- 
tween the electron and photon equations can be simplified because the higher cascades have 
negligible influence in this energy region, and (3) it is possible, though not necessarily desirable, 
to seek analytic solutions. 

Because even an analytic development leads to a numerical integration, it appeared most 
expeditious to employ eq (3), with the simplifications just mentioned and a finer integration 
spacing, to develop a solution in this region. In particular, for an electron-generated cascade, 
the term for generating electrons from photons and the term A[dF(E)/dE] are very small and 
can be neglected. The remaining functions K can be written in the approximate form: 
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K°°(E',E)=z^-B\n 



E'-E 

E 



■0, 



tf«{E) =A In ^^-B(E ~E) In E "~ E 



Q 

where A, B, and C are constants (see appendix 1). 



eE 



■C(E Q -E), 



The quantities (7) were tabulated at E X2 w/6 , where m=^ -> ^> 1. 

o 4 Z 



This spacing required 



a special set of weights; but otherwise the evaluation in this special case is carried out by 
using the appropriate form of eq (5) in a manner similar to the main calculation. Results 
of this special solution not only give the solution near E but also permit a more accurate 
evaluation of the integrals in eq (3) for a,E f >E X2~ l/Q . In this energy region the use of 
Simpson's rule with 2~ 1/6 spacing is not really a good enough approximation because F(E) 
has large nonquadratic variations. 

8. Calculations 

The calculations in this paper are based on a monoenergetic source of 1 electron per second 
at 360 mc 2 in lead. Results were first obtained for the special solution near E . The main 
calculation, extending down to 5.6 mc 2 , was carried out in two parts, using the mesh E X2~ m/& 
to the energy \E Q and the mesh E X2~ m/ * from that energy to E=5.6 mc 2 . 





Figure 1. Arrangement of tabulated functions. 



The calculation was facilitated by arranging a triangular tabulation of the four W's 
and a linear tabulation of the C's and Dina form shown in figure 1 . Most of the work went 
into tabulating K functions and W's, the final computation being carried out in a few hours. 
The results are given in table 1 . 

Any error had the effect of introducing an oscillation in the final computation, with alternate 

points being too high and too low. The integration procedure appeared stable in that such 

oscillations tended to damp out slowly rather than to magnify. By plotting the results it was 

possible to use these oscillations to find errors. The final sequence of numbers is smooth, 

i. e., nonoscillatory, to better than 1 percent, which may indicate that the integration procedure 

introduced errors no larger than this, although the possibility that errors may have piled up 

at low energy is not excluded. 

Table 1. Results 



E 


F(E) 


H{a) 


E 


F(E) 


H(«) 


E 


F(E) 


H(a) 


360.0 
354.8 
349. 8 
339.8 
320. 7 

285.7 
254.6 
226.8 
202.0 
180.0 




0.751 
.980 
1.35 
1.76 

2.47 
3.05 
3.61 
4.14 
4.71 





160.4 
142.9 
127. 3 
113.4 
101.0 

90.0 
71.4 
56.7 
45.0 
35.7 


5.26 
5.92 
6.52 
7.31 

7.97 

8.91 
10.5 
13.0 
15.6 
18.9 


0. 0409 
.0537 
. 0693 
.0901 
.115 

.148 
.235 
.384 
.601 
.944 


28.4 
22.5 
17.9 
14.2 
11.3 

8.93 
7.09 
5.63 


22 A 
26.8 
31.4 
36.5 
41.6 

47.9 
53.8 
60.3 


1.44 
2.20 
3.27 
4.80 
6.85 

10.1 
13.8 
19.2 




0. 00488 
. 00652 

. 00985 

.0135 

.0180 

.0238 

.0314 
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9. Discussion of Results 

The results in table 1 are plotted in figures 2 and 3, where they are compared with the 
results of Richards and Nordheim [6] and calculations based on the analytic solution given 
by Rossi and Greisen [7]. (See appendix 3 for a discussion of the normalization.) For the 
comparatively low initial energy considered here, the equations of Rossi and Greisen are only 
accurate in a very limited energy region around 130 mc 2 . At lower energies, as pointed out 
by Richards and Nordheim [6], a more accurate representation of the interaction probabilities 
is necessary. The Richards and Nordheim numerical calculation was based on more accurate 
cross-sectional forms, but their results are normalized to those of Rossi and Greisen at 131 mc 2 . 
Thus, near that energy, the significant comparison is with Rossi and Greisen. At lower 
energies, due to differences in the cross-sectional functions involved, the three curves go aparl . 
At energies near E no comparable calculations have been found. Thus the only energy 
region where close agreement between the spectra calculated in this paper and those given by 
the other authors could be expected is from % to \E$. In this region the agreement is shown 
in figures 2 and 3 to be quite good. 

Through a consideration of the specific cross-sectional forms used in the three calcula- 
tions, some further comparisons are possible at lower energies. For example, the spectra of 
Rossi and Greisen indicate a larger number of electrons and a smaller number of photons than 
reported in this paper. This is caused mainly by the cross section for pair production, which 
is significantly lower under the complete screening approximation than the asymptotic form 
used by Rossi and Greisen. Thus, in the calculations in this paper, low-energy photons do 
not generate as many electrons, and there results a decrease in F{E) and an increase in 11(a) 
relative to Rossi-Greisen. The spectra of Richards and Nordheim indicate a much larger 
number of electrons and photons at low energies. Although this must be traceable similarly 
to a difference in cross sections, the reasoning is not as straightforward in this case. 
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Figure 2. Comparison of F(E) with previous results. Figure 3. Comparison of H (a) with previous results. 
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10. Final Comments 

The agreement as to general trends and as to absolute magnitude near %Eq demonstrates 
the feasibility of the calculational procedure developed in this paper. Logically, the next 
step is to introduce accurate cross-sectional forms into the equations to obtain more accurate 
spectra. Alternatively, one could move immediately toward a solution of the spatial and 
directional distributions. This would involve computing the higher moments in the linked 
set of equations of which these spectra form the zeroth moment (for an example of this, see 
reference [10]). 

11. Appendixes 

Appendix 1 . Cross-Sectional Forms 

The following detailed cross-sectional forms were used to derive the functions K discussed 
in section 5. 

Bremsstrahlung losses [11]: 

„<E',E> -45 ± (ji) [(i±i-|) In 18 3Z-» + i} 

where rj=E/E\ 

Electron atomic interactions [8] : 

cr M {E',E)=\ <p T -^ [^-Lp+I], 

where 7]=E/E'. 

Pair production (either electron or positron) [11]: 



S,E)=8f ± [(| „ 2 -| ,+l)ln 183Z-»-i *(l-l)]» 



where rj = E/a\ 

Compton electrons: 

Compton photons: 



( , , 3 i rv+n 



where rj=a/a f . 

Bremsstrahlung photons [11]: 



c bp {E',a)=^>\ [(|-| V + V ^ In 183Z-H+i^f], 



where rj=a/E\ 

The cross-sectional constants were evaluated as follows: 



and 



v= 137 ~^= 3 - 9435 barns > Z=82; 



<Pt=o t?"o= 0.6652 barn. 



The constants A, B, and C in the special solution near E (section 7) are obtained from the 
expression for K es (E',E); they are coefficients of the terms proportional to E~ l (l — 77)" 1 , 
— ln(l— 7}), and 1, respectively, in the limit E,E'-^E : 
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.4=A T Z% r =0.1189 



i?=A74^f"^ln 183Z-^+i"|=0.2339 
C=mj [-£ In 183Z-*-i]-f NZ, T |=- 



-0.1487. 



Appendix 2. Integration and Differentiation Weights 

To obtain the integration weights u in used in conneetion with eqs (5) and (6), the following 
formulas were used : 

Jx -A 6 

( X{&A dxf(.r) « A [5/fo + A) +8/(z ) -y^o - A) ] , 
f^Wto «^ [f(x +2A)+Zf(x Q )], 

J.r -A 4 

where /(a) is a function that varies slowly for x in the region of integration. 

For differentiation the weights a t were determined by using tin 4 following formula: 



0/M„) 
d/'J 






■u(E ll -E J y 



)*i 






Appendix 3. Normalization of Previous Results 

The following results were obtained by using a table of calculations given by Richards 
and Nordheim [6]. 



E 


F(E)=mMF RN {E) 


H(a)=5.S29g(K) 


131 


6. 43 


o. 0639 


91. 8 


9. 08 


. 171 


65. 6 


12. 6 


. 384 


39. 3 


21. 


1. 12 


26. 2 


30. 


2. 61 


19. 7 


37. 7 


4. 74 


13. 1 


49. 6 


10. 7 


6. 56 


71. 3 


34. 1 



The normalization factors arise in the following manner: 



69.86=|0.437X~^lxdensity of lead (11. 005) Xradiation length (0.529) 



5.329 = 69.86 X energy conversion factor 



13.11 



For a discussion of the factor in square brackets, see reference [6]. 
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The following results were obtained by using the equations of Rossi and Greisen [7] under 
"approximation B": 



[10] 
[11] 



E 


F(E) = b.822Z T (E ,E) 


H(a)=2.975z y (Eo,W) 


180 
90 
45 
22. 5 


4. 77 
9. 02 
16. 4 
28.5 


0. 0348 
. 132 

. 482 

1. 70 



The normalization factors arise in the following manner: 

5.822 = density of lead (11.005) Xradiation length (0.529), 
2.975=5.822Xmc 2 (0.511). 
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